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All-dielectric ” magnetic light” nanophotonics based on high refractive index nanoparticles allows 
controlling magnetic component of light at nanoscale without having high dissipative losses. The 
artificial magnetic optical response of such nanoparticles originates from circular displacement cur¬ 
rents excited inside those structures and strongly depends on geometry and dispersion of optical 
materials. Here a new approach for increasing magnetic response via resonant bianisotropy effect 
is proposed and analyzed. The key mechanism of enhancement is based on electric-magnetic in¬ 
teraction between two electrically and magnetically resonant nanoparticles of all-dielectric dimer 
nanoantenna. It was shown that proper geometrical arrangement of the dimer in respect to the 
incident illumination direction allows flexible control over ah vectorial components of magnetic po¬ 
larizability, tailoring the later in the dynamical range of 100% and enhancement up to 36% relative 
to performances of standalone spherical particles. The proposed approach provides pathways for 
designs of all-dielectric metamaterials and metasurfaces with strong magnetic response. 


INTRODUCTION 


Intrinsic magnetic polarizabilities of natural materi¬ 
als have strong frequency dependence with the funda¬ 
mental cut-off in GHz range, originating from relatively 
low spin and orbital susceptibilities [1]. Recently, effec¬ 
tive polarization currents in subwavelength structured 
loops, organized in ordered arrays, became sources of 
high-frequency artificial magnetism [2]. Nanostructured 
noble metals, supporting localized plasmon resonances 
could serve as building blocks for metamaterials with 
artificial magnetic polarizability [3]. However, inherent 
material losses set severe limitations on performances of 
such structures Q. Another approach for obtaining mag¬ 
netic optical response is to employ circular displacement 
currents in high-index dielectric nanoparticles [5j . This is 
the essence of so-called all-dielectric nanophotonics which 
opened the way to control magnetic component of light 
at nanoscale without high-dissipation, inherent for metal¬ 
lic nanostructures [6H12|. The ” magnetic light” concept 
found use in various applications, such as nanoanten¬ 
nas 00, quantum interface for NV-centers 0 , pho¬ 
tonic topological insulators M > broadband perfect re¬ 
flectors 171 7 w aveguides ||18l], cloacking 0 , harmonics 
generation [20|, wave-front engineering, and dispersion 
control [211. 

Magnetic response of a dielectric particle strongly de¬ 
pends on it’s refractive index, shape, and external envi¬ 
ronment. The eigen frequencies of electric and magnetic 
resonances could span the entire visible range. However, 
the value of those multipole moments is limited by dis¬ 
persion properties of optical materials. 


Here a new approach for tailoring magnetic response of 
dielectric nanoparticles via resonant effect of bianisotropy 
is proposed. Microscopically, bianisotropy is the effect 


of magneto-electric coupling, where electric polarization 
induces magnetic and vice versa. The signature of this 
effect appears in the constitutive relations, e.g the depen¬ 
dence of electrical induction also on magnetic field and 
magnetic induction also on electric field 22|, [23j. Bian¬ 
isotropy is used for achieving high values of effective po¬ 


larizability in metamaterials [24j], unique properties of 
metasurfaces |23|, and designed directivity of nanoanten¬ 
nas [ii] . Previously, it was shown that the interaction of 
dielectric nanoparticles with substrates may increase in¬ 
duced magnetic moment due to the effect of non-resonant 



FIG. 1. The geometry of the bianisotropic all-dielectric 
dimer nanoantenna. The nanoantenna consists of two dielec¬ 
tric nanoparticles separated by a distance r 2 i. The sizes of 
nanoparticles are Ri (bigger nanoparticle) and R 2 (smaller 
nanoparticle). The bigger nanoparticle exhibit electric dipole 
resonance, while the smaller one magnetic dipole resonance at 
the same wavelength (480 nm). The nanoantenna is excited 
by a y-polarized plane wave propagating along a>axis. 
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bianisotropy (2(|. The approach proposed here is based 
on electric-magnetic interaction between two resonant 
nanoparticles of dimer nanoantenna (see Fig. [1]). The 
nanoparticles were designed in the following way: the 
electric dipolar resonance of the bigger sphere overlaps 
with magnetic response of the smaller one. In this case, 
the effect of bi-resonance anisotropy is achieved: the res¬ 
onant electric moment of the bigger nanoparticle induces 
the additional magnetic moment in the smaller one, tai¬ 
loring its overall response. 

The manuscript is organized as follows: first, the op¬ 
tical properties of isolated spherical particles are briefly 
discussed in the context of the resonance tuning. Next, 
the analytical coupled dipoles formulation of the problem 
is developed and verified by numerical modeling. The ex¬ 
pression describing the magnetic moment of the nanopar¬ 
ticle, considering the bianisotropy, effect is derived. Fur¬ 
thermore, it was shown that proper geometrical arrange¬ 
ment of the dimer in respect to the incident illumination 
direction allows achieving additional vectorial component 
of magnetic polarization. 

COUPLED DIPOLES THEORY AND 
NUMERICAL RESULTS 

In order to obtain the dimer design, properties of iso¬ 
lated components will be briefly discussed. First, iso¬ 
lated silicon sphere of radius R 2 = 52 nm, having mag¬ 
netic dipolar resonance at wavelength 477 nm, i.e. in the 
visible range is considered. The material dispersion of 
crystalline silicone (c-Si) is taken from [ 27 }. The scat¬ 
tering cross-section of the nanoparticle have been calcu¬ 
lated using CST Microwave Studio, and corresponding 
results as the function of wavelength are presented in 
Fig. (2fa) (blue curve). Electric dipolar resonance is blue 
shifted in respect to the magnetic one and appears on the 
spectrum at 418 nm wavelength. This Mie resonances 
hierarchy is red shifted when the radius of nanoparti¬ 
cle is increased. It is worth noting, that high-order 
quadrupole resonances, also contributing to the scatter¬ 
ing cross-section, are suppressed in this region due to per¬ 
ceptible losses of silicon. The bigger sphere with radius 
of R\ = 70 nm exhibits its electric resonance at the sim¬ 
ilar spectral position where 52-nm sphere exhibits reso¬ 
nant magnetic response (Fig. [2^, red curve). The electric 
and magnetic polarizabilities which are associated with 
induced electric and magnetic moments have been also 
calculated (Fig. (2)3). Results of numerical simulations 
(scatter plots) are in good agreement with analytical Mie 
theory (lines plot) verifying the validity of the numerical 
tool. Values of those electric and magnetic moments will 
be subsequently used in the analytical model, based on 
discrete dipoles approximation. 

Next, the scattering of a plane wave on all-dielectric 
dimer nanoantenna will be analyzed. The electromag- 
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FIG. 2. Optical properties of single silicon nanoparticles with 
radii R\ — 70 nm (red curves) and R 2 = 52 nm (blue curves), 
(a) The scattering cross-section spectra, normalized to geo¬ 
metric cross-section (nr 2 ), (b) Dispersion of particles polar¬ 
izabilities. The curves correspond to analytical calculations, 
red circles and blue squares represent numerical results. Elec¬ 
tric dipolar resonance of 70-nm radius particle overlaps with 
the magnetic dipolar resonance of 52-nm sphere at A ~ 480 
nm. 


netic scattering problem could be solved by employ¬ 
ing Coupled Electric and Magnetic Dipole Approxima¬ 
tion [28[ (CEMDA). In this method, complex nanostruc¬ 
tures are represented by converging series of point elec¬ 
tric and magnetic dipoles [29], while the problem of two 
spheres could be approximated by only two. This ap¬ 
proximation is particularly accurate, if the gap between 
the spheres is bigger than their radii [30|. The goal of 
the subsequent analytical analysis is obtaining a simple 
formulation, underlining the interference phenomena, af¬ 
fecting the magnetic dipolar polarizability of the smaller 
particle and revealing the bianisotropic nature of the in¬ 
teraction. The high-order modes and their interplay have 
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been taken into account using the full-wave numerical 
calculations using the CST Microwave Studio. Following 
the CEMDA method, total electromagnetic fields [both 
electric (E) and magnetic (H)] were decomposed into in¬ 
cident and scattered components: 


E(r,) = E 0 (r i ) + Ef(r,), 
H(r i )=H 0 (r i )+Hf(r i ), (1) 

where indices = 1,2 (i % j ) denote first (bigger) 
and second (smaller) nanoparticles, respectively, Eo and 
H 0 are electric and magnetic fields of the incident plane 
wave, E(ri) = E(0) and H(ri) = H(0) are the full elec¬ 
tric and magnetic fields at the position of the first (elec¬ 
tric) nanoparticle (with radius-vector 1*1 = 0), E| c (ri) 
and H^iq) are the electric and magnetic fields scat¬ 
tered by the second (smaller) nanoparticle at the first 
particle center, and analogously for the second particle. 
Both particles are polarized by the incident field, as well 
as by the scattered one: 

Pi = £ 0 af E(r i), m; = a 8 M H(rj), (2) 


where af and af 1 are the electric and magnetic polariz¬ 
abilities of single particle ” i”: 


E . 67 T£h m • 

a i a i =lT3 bl ’ 

K h K h 


( 3 ) 


and eo is the dielectric constant, j is the imaginary 
unit. The coefficients a\ and b\ are called the first or¬ 
der Mie scattering amplitudes, and they can be writ¬ 
ten in simplified form which is suitable for analysis: 
ai(X) = [. A-B]/[C-D ], 61 (A) = [ Bn- 2 -A\/[Dn~ 2 -C ], 
where the following notation is introduced: 
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where n = y / e/sh , p = khRi, R% - radius of the sphere 
”i”, e - the permittivity of the material (silicon, here), 
kh = 2iT^Sh/\ is the wave number in the surrounding 
space, and A is a wavelength. 

The values of single particle polarizabilities, calculated 
using Eqs. m and are in a perfect agreement with 
our numerical calculations [see Figl2][b)]. Thus, the scat¬ 
tered fields of the dipoles can be obtained through the 
Green function of a point dipole in free-space G(r^, iq) = 
Gij UK, following the method CEMDA: 


Pi = a* 


E 0 (i*i) + — GijPj + -y—[g ij x m 


^0 


m i = a 
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ck 0 
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where fco and c are the wavenumber and speed of light 
in vacuum. The Green’s function of a point dipole in 
free-space is well known: 
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where r^ — Yj is the radius vector, connecting the 

center of the first dipole (coordinate origin, or center of 
the bigger particle) with the second one (center of the 
smaller particle). The tensor and vector have the 
following symmetry of indices permutation: Gij = Gji 
and gij = -gji. 

Both nanoparticles, being isolated, have 3-fold degen¬ 
erated (magnetic and electric) dipolar resonances, ori¬ 
ented along the unit vectors of a Cartesian coordinate 
system. The excitation of an isolated sphere is solely 
defined by the polarization of the incident wave - for ex¬ 
ample, linearly polarized beam will excite only one of 3 
components of the dipolar mode. However, the geometry 
of the coupled dipoles together with the excitation (not 
necessarily coinciding with one of the symmetry axis of 
the system) will break the degeneracy and, as the re¬ 
sult, all three vectorial components must be taken into 
, account. The obtained set of equations can be solved an¬ 
alytically by means of the matrix inversion or numerically 
\ in the same fashion. In order to verify the validity of the 
) ’ proposed theoretical model, we consider a particular case, 
where the system of Eqs. ([5]) has a simple and intuitive so¬ 
lution. Arranging the nanoparticles along the x-axis and 
exciting the system with linearly polarized plane wave 
( 4 ) along y-axis with angles being 0 = 7 r /2 and tp = 0 , the 
set of 12 coupled equations was reduced just for 4, since 
the symmetry considerations allow the moment compo¬ 
nents to be induced only along z-axis (magnetic dipole) 
and y-axis (electric dipole). Similar configuration was 
studied in [31( for hybrid metal-dielectric nanoantennas 
with non-resonant dielectric nanoparticles. This set of 
coupled equations has particularly simple solution, which 
agrees well with the numerical calculations. The values 
of the magnetic polarizability enhancement, calculated 
both analytically and numerically as the function of the 
distance 1*12 between nanoparticles are shown in Fig. 02 
It is clearly seen, that the analytical and numerical mod¬ 
els agree with each other and showing 15% enhancement 
of the magnetic moment of the single nanoparticle, as 
the result of the bianisotropic coupling. For clear un¬ 
derstanding of oscillatory behaviour of the magnetic mo- 
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FIG. 3. Enhancement of magnetic polarizability of the smaller 
nanoparticle for the case 0 = n/2, p = 0. Red curve corre¬ 
sponds to the solution Eq. ©, black squares are the results 
of numerical full-wave simulation. 

merits, the simplest model, when the polarizabilities af 1 
and a f are supposed to be zero, was considered. In this 
case the simple formula can be obtained in form showing 
the bianisotropic nature of the effect: 

m 2 = t]Hq + 7 ^ 0 , 

a M e jk 0 r 21 

^ ~ 1 _ rv E rv M n 2 P’ 

i aq a 2 gijK 0 

_ -eoafa^gijjkoc 

7 " l-afo «g%kl 

The qualitative analysis of magnetic moment m 2 enables 
to observe the clear interference phenomena - the direct 
excitation of magnetic dipole by the plane wave [alf e lk ° r 
term in Eq. ©] and the contribution of the scattered field 
through the term a 2 a^ . Thus, oscillation behavior of 
magnetic moment enhancement is the result of the inter¬ 
ference phenomenon (see Fig. [3]). For longer distances 
the coupling between the particles becomes weaker and 
converges to the value of the isolated particle. An almost 
perfect fit of full numerical modeling with the CEMDA 
method enables to use the later for analysis of more com¬ 
plex structures without involving full wave simulations. 
This approach will be subsequently employed. 

VECTORIAL STRUCTURE OF MAGNETIC 
MOMENTS 

In the subsequent studies the illumination was cho¬ 
sen to propagate along x-axis and being polarized along 
y. There are three geometrical parameters, affecting the 
magnetic polarizability of the smaller particle: the dis¬ 
tance between the spheres’ centres r 2 \ = Ri+R 2 +D, and 


D=100nm D=200 nm D=300 nm 
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FIG. 4. The three-dimensional angular dependencies of the 
relative increase of the components l = x,y, z of the smaller 
nanoparticle magnetic moment ({m^-^/a^ 1 Ho\) for the differ¬ 
ent distances D=100, 200, and 300 nm. 


the angles 9 and ip. The angular dependence of the in¬ 
duced magnetic moment components of the smaller parti¬ 
cle for 3 characteristic separation distances D = 100, 200, 
and 300 nm (gap between the nanoparticles’ surfaces) will 
be studied next. In the case of noninteracting nanoparti¬ 
cles there is the only one non-zero component of nanopar¬ 
ticle’s magnetic moment - codirected with the mag¬ 
netic field. In the case of bianisotropic coupling there are 
also x- and y- non-zero components of the magnetic mo¬ 
ment (m^, m^). Fig. |4] shows the three-dimensional 
angular dependencies of normalized magnetic moment 
components (\m 2 '^l = x,y,z) of the smaller 
nanoparticle, where alf is the polarizability of the sin¬ 
gle smaller nanoparticle. These results have been ob¬ 
tained by exact solving Eqs. 0 taking into account 
the electric and magnetic responses of both nanoparti¬ 
cles. As it could be seen, the induced components are 
smaller than the main moment. Those additional vec¬ 
torial components of magnetic moments increase when 
the distance between the particles goes down and may 
being up to 50% for D < 50 nm. But it should be 
noted that for such small distances between the parti¬ 
cles accuracy of CEMDA method is low, sow this value 
is rough [30[. When the distance between the particles 
goes up, secondary vectorial components of magnetic mo¬ 
ment decrease. Therefore, the overall variation of mag¬ 
netic moments is mainly determined by its ^-component. 
Diagrams of ^-component of the magnetic moment are 
asymmetrical and their forms are dissimilar for differ¬ 
ent separation distances between the particles (D). For 
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Enhancement of magnetic moment, % 



Wavelength, nm Wavelength, nm Wavelength, nm 

FIG. 5. The spectral dependencies of the enhancement of magnetic polarizability [{jm 2 /a^ Ho\ — l) • 100%] of the smaller 
nanoparticle for D — 168 nm and for different angles (tp,0) on wavelength. The data for 0 G (tt/2, 7r) and ip G ( 7 r, 27 r) repeats 
the presented results, therefore they are omitted. For the case of 0 = 0 the lines for all ip are coincide. 


D > 300 nm enhancement of magnetic moment is weak, 
because of the vanishing coupling between the particles. 
For D >2 fi m this diagram is symmetrical with good ac¬ 
curacy, as it replicates the performance of the isolated 
particle. For smaller distances D the maximum of mag¬ 
netic moment corresponds to the cases (0 = 7t/2, cp = 7t/2 
and ip = 37t/2), and for bigger distances maximum of 
magnetic moment corresponds to the case (0 = 7t/2, 
(p = 7r). It should be noted that case (0 = 7t/2, (p = it/2 
and ip = 37t/2) corresponds to magnetic dipoles coupling 
only, while the case (0 = 7t/2, ip = 7 r) corresponds to 
electric-magnetic dipoles coupling only. Additional anal¬ 
ysis shows that the maximum of relative magnetic mo¬ 
ment enhancement is nearly 36% and it is achieved for 
D = 168 nm for the case of (0 = 7r/2, cp = 7 r). The 
performance of the dimer with this separation distance 
will be investigated in details hereafter. Both induced 
magnetic and electric moments of the smaller nanoparti¬ 
cle were calculated using the theoretical model CEMDA. 
Enhancement of the magnetic polarizability of smaller 
nanoparticle (| mn^l H§\ — l) • 100% for different angu¬ 
lar arrangements (the distance D =168 nm is kept con¬ 
stant) is shown in the Fig. [5j These dependences show 
that the effect has a resonance character. Moreover, the 
maximum value is achieved at the wavelength of 480 nm 
i.e. at the electric dipole resonance of bigger nanoparticle 
and magnetic dipole resonance of smaller one underlining 
the impact of the resonant nature. 

OUTLOOK AND CONCLUSIONS 

Coupled particles approach for controlling magnetic 
polarizabilities of nanoscale spheres was proposed. While 
standalone nanoparticles allow obtaining dipolar and 
high-multipolar resonant responses at the desired wave¬ 
length in the visible range specifying their radii and ma¬ 
terials, a system of two coupled nanoparticles possesses 


more degrees of freedom. Altering the radii of both 
nanoparticles allows to investigate the impact of all the 
combinations of multipolar coupling effects on the prop¬ 
erties of the system. Naturally, the amplitude of the 
effects is dependent on their mutual displacement, and it 
decreases at longer distances between the particles, as it 
was proved both analytically and numerically in the case 
of dipolar magnetic-electric coupling. Symmetry consid¬ 
erations allow exciting only one dominant induced elec¬ 
tric and magnetic moment component when a plane wave 
is incident upon a single nanoparticle. It was shown, that 
for a nanoparticle system, altering the spherical angles 
and distance between nanoparticles allows to excite all 
vectorial components of moments simultaneously. The 
secondary components being up to 50% of the amplitude 
value of dominant ones. Furthermore, these parameters 
also define the spectral position of deeps and peaks of 
the dominant electric moment components and ampli¬ 
tude values, allowing full on-demand control of the elec¬ 
tromagnetic properties of the system of coupled nanopar¬ 
ticles. The proposed approach can find use in designs of 
more complex structures such as all-dielectric metamate¬ 
rials and metasurfaces with strong magnetic responses. 
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